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We introduce a group naturally acting on aperiodic necklaces of length n 
with two colours using the 1-1 correspondences between aperiodic necklaces and 
irreducible polynomials over the field F2 of two elements. We notice that this 
group is isomorphic to the quotient group of non-degenerate circulant matrices 
of size n over that field modulo a natural cyclic subgroup. Our groups turn 
i-Q . out to be isomorphic to the sandpile groups for a special sequence of directed 

^ \ graphs. 
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1 Introduction 

This work originated in the research of the first author related to the Drinfeld asso- 
10 ! ciator |DD] . It is well-known (see, e.g. |CDM] ) that the logarithm of the classical 

^^ ' associator ^kz is an element of the completed free Lie algebra on two generators 

with coefficients in the algebra Z of multiple zeta values [Hoff j . The free Lie algebra 
Cf) • has a basis whose elements are labelled by aperiodic necklaces |Reu] . However, the 

explicit expansion of log ^kz over this basis shown in [DD] displays a highly chaotic 

behaviour. Therefore, a natural idea arises to introduce some structure in the set of 
^ I necklaces. In this paper, we introduce a group acting on the set of aperiodic necklaces 

of fixed length in the hope that the orbits of this group may shed some light on the 

structure of the embarrassing expression given in |DD] . 

In the spirit proclaimed by V. I. Arnold in |Arn] , this paper does not contain any 

proofs, only constructions, problems, motivations, examples and statement of results 

proven elsewhere. It may be considered as an informal introduction to the paper 

PlP] . 

2 General setting 

In the most general setting, the approach we follow here consists in the following. 
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Suppose we have two sets A and B and a family of bijections ipi : A ^ B. This 
data can be used to define two groups: Ga, acting on the set A, and Gb, acting on 
the set B. Indeed, let Ga be the group generated by all bijections ipij = ip~^ o ip^ and 
Gb the group generated by all Lp[^ = fi^ V^J^- 

The following lemma is immediate. 

Lemma 1. The groups Ga and Gb are isomorphic, an isomorphism being given by 
the assignment g \-^ (pi o g o ip~-^ for any fixed i. In particular, the actions of these 
groups on the sets A and B are equivalent, any bijection cpi maps the orbits of the 
group Ga onto the orbits of the group Gb, and this map on the set of orbits does not 
depend on the choice of a particular i. 

3 History 

The story began when the first author obtained the results of Section |6] (and also 
Section [TT] which is, however, irrelevant to the main topic of this paper). This was 
done experimentally, by computer. The second author found the sequence of orders of 
the groups RG^ in Sloane's encyclopaedia (A027362) and conjectured that the groups 
are isomorphic to the sandpile groups of generalized de Bruijn graphs. This conjecture 
was checked by computer up to order 16 by the present authors and then proved by 
S.H.Chan in his Bachelor's thesis |Ch] IH Finally, this theorem was generalized and 
extended in the paper jCHPj . 

4 Groups acting on necklaces 

A p-coloured necklace of length n is a sequence of n objects of p different kinds (called 
beads) considered up to cyclic shifts. Necklaces may be periodic (admitting a non- 
trivial shift that does not change it) or aperiodic. If the colours are linearly ordered, 
then the lexicographically minimal of all cyclic shifts of an aperiodic necklace is called 
a Lyndon word. We will always suppose that p is prime and use the elements of 
a prime finite field ¥p as the colours of beads. We will refer to such necklaces and 
Lyndon words as arithmetical. 

The main case which is most interesting from the point of view of Section [T] is 
p = 2. 

Example: There are exactly 6 aperiodic 2-coloured necklaces of length 5 given by 
Lyndon words 00001, 00011, 00111, 01111, 00101, 01011. 

Denote the set of aperiodic necklaces with parameters n and p by A^^. There is a 
formula for the cardinality of this set in terms of the Mobius function: 



l^nl = ^EM^)P" 



n/d 

n 



d\ri 



^Actually, the correct proof appeared only in the updated version of the thesis paper. 



It is remarkable that this number is equal to the number of irreducible polynomials 
over the field ¥p of degree n, or, which is the same, to the number of orbits of maximal 
length n of the action of the Galois group Gal(Fpn/Fp). □ Denote the set of irreducible 
polynomials by /^. There are two known explicit constructions of a 1-1 correspondence 

T ■ ^ ' n n 



One of them (first mentioned in Reutenauer |Reu] who ascribed it to folklore) depends 
on the choice of a normal polynomial of degree n over Fp. Another one belongs to 
Golomb |Gol] and depends on the choice of a primitive polynomial of degree n over Fp. 
According to Section Ej either set of bijections generates a group of transformations 
on the set of aperiodic necklaces. We will call the first one the Reutenauer group and 
denote it by -RG^, and the second, the Golomb group and denote it by GG^. 



5 Galois group and irreducible polynomials 

It is well known that the Galois group Gal(Fpn : Fp) is a cyclic group of order n gen- 
erated by the Frobenius automorphism o"(x) = x^. Each orbit of this group coincides 
with the set of roots of an irreducible polynomial over Fp of degree which divides 
n. The orbits of maximum length n correspond to irreducible polynomials of degree 
exactly n. The union of all such orbits is equal to the complement in Fpn of all its 
proper subfields. A polynomial is called normal, if the set of its roots constitutes a 
basis of the vector space Fpn over Fp; it is called primitive, if one (and hence any) of 
its roots is a generator of the multiplicative group F*„. 

Here is an example. Consider the extension Fie '■ ^2 of degree 4 as the quotient 
ring F2[X]/(X^ + X + 1). Denoting the class of X in this quotient by a, we can see 
that a is a generator of Fie, and the following table lists the orbits of the Galois group 
Gal(Fi6 : F2), the corresponding irreducible polynomials and, for the polynomials of 
maximum degree, indicates their nature: 



orbit 


polynomial 


normal? 


primitive? 


{0} 


X 






{1} 


x + 1 






K,«^°} 


x'^ + X + 1 






{a,a'^,a'^,a^} 


x'^ + x + 1 


no 


yes 


{a\a^a'^,a'} 


x^ + x^ + x"^ + x + 1 


yes 


no 


{a',a'\a'\a''} 


x^ + x-^ + l 


yes 


yes 



^This is why we consider only the necklaces with a prime number of colors. Generalizing the 
theory for p'' instead of p is an interesting open problem. 



6 Reutenauer's construction 

By definition, a normal basis is an orbit of tlie Galois group which constitutes a basis 
of Fpn as a vector space over ¥p. Given a normal basis A = {a, a^, a^ , • • • } and a 
necklace z/q, z/i, . . . , fn-ii we set up a sum z/qO + z/io^ + ■ — h Un-iO^"^ . Cyclical shifts 
on the sequence z/q, z/i, . . . lead to the change of the resulting element of the big field 
Fpn within the same orbit of the Galois group, so that the minimal polynomial of that 
element remains the same. 

We have therefore two finite sets of equal cardinality, N^ and /^, equipped with 
a family of 1-1 maps ipA '■ N^ — )■ I^. It turns out that the group generated by this 
set according to the construction of Section |2l coincides in this case simply with the 
set of all maps (f^^ o (^^ of the set of necklaces into itself. It forms an abelian group, 
earlier denoted by RG^, whose order is equal to the number of normal bases (in the 
case p = 2 known as Sloane's sequence A027362). 

Here is the table of automorphism groups of necklaces for p = 2 and n < 15 (M„ 
stands for the number of normal bases, that is, the order of the group, and the last 
column lists the lengths of the orbits of RG^ on the corresponding set of necklaces) : 



n 


m 


M„ 


group 


orbits 


2 


I 


1 


Zi 


1 


3 


2 


1 


Zi 


1 + 1 


4 


3 


2 


Z2 


2 + 1 


5 


6 


3 


Z3 


3 + 3 


6 


9 


4 


Z2 0Z2 


4+2+2+1 


7 


18 


7 


Z7 


7+7+1+1+1+1 


8 


30 


16 


Z2 © Z2 © Z4 


16 + 8 + 4 + 2 


9 


56 


21 


Z21 


21 + 21 + 7 + 7 


10 


99 


48 


Z2 © Z2 © Z2 © Ze 


48 + 24 + 24 + 3 


11 


186 


93 


Z93 


93 + 93 


12 


335 


128 


Z2 © Z2 © Z2 © Z4 © Z4 


128 + 64 + 32 + 32 + 16 + 16 

+ 16 + 8 + 8 + 4 + 4 + 4 + 2 + 1 


13 


630 


315 


Z315 


315 + 315 


14 


1161 


448 


Z2 © Z2 © Z2 © Z2 © Z2 

©Zi4 


448 + 224 + 224 + 56 + 56 

+28 + 28 + 28 + 28 + 8 + 8 + 7 
+4+4+4+4+1+1 


15 


2182 


675 


Z3 © Zi5 © Zi5 


675 + 675 + 225 + 225 

+45 + 45 + 45 + 45 + 45 + 45 
+ 15 + 15 + 15 + 15 + 15 + 15 
+3+3+3+3+3+3 
+1+1+1+1 



7 Orbits of the Reutenauer group 

It is very interesting to study the orbits of the group RG^ acting on the set of aperiodic 
necklaces N^. So far, we only have some empirical results in the case p = 2. 
For example, 

• If n = 4, the orbits are Oi = {0001, 0111} and O2 = {0011}. The action of the 
group -RG4 = Z2 on each orbit is evident. 

• If n = 5, we have Oi = {00001, 00111, 01011} and O2 = {00011, 01111, 00101}. 
The action of the group is cyclic on each orbit. 

• If n = 6, then d = {000001,011111,001011,001101}, O2 = {000011,010111}, 
O3 = {000101,001111} and O4 = {000111}. Here, the automorphism group 
acts as Z2 © Z2 on Oi, as Z2 on O2 and O3, and trivially on the last orbit. 

Returning to ideas of Section [H we tried to evaluate various symmetric functions 
of the coefficients in the Drinfeld associator over the orbits in these cases, but could 
not arrive to any sensible conjecture. 

It is also worthwhile no notice that there is an interesting operator on the set of 
necklaces, which sometimes turns an aperiodic necklace into a periodic one, but in 
general it takes the whole orbits in the above lists into whole orbits. We call it the 
averaging operator; by definition, it acts as follows: {z/j} 1— )■ {vi + Vi+i}- For the above 
examples, we have 

• n = A: Oi — VO2 — ^ 0. 

. n = 5: Oi — y02 — > O2. 

• n = 6: Oi ^ O2 ^ O3 ^ O3, O4 — ^ 0. 

(Here, going to means that the necklace becomes periodic.) 

We noticed that in all examples with n < 15, there is a main orbit, that is, an 
orbit of maximal length equal to the order of the group RG^ which acts on this orbit 
simply transitively. Iterated averaging operators applied to the main orbit give the 
majority of orbits; some smallest orbits may go to 0. 

Remark that here we spoke about the orbits of arithmetical groups in the sense 
of Section |H From the point of view of studying the Drinfeld associator, however, it 
makes little sense to distinguish between symbols and 1: it is more reasonable to 
extend our groups RG^ by the operator that fiips O's and I's in each necklace. The 
orbits of such extended groups consist of either one or two orbits of the initial groups; 
it should be interesting to have their explicit description and try to study, for example, 
the sums of coefficients in the logarithm of the associator over these extended orbits. 



8 Sandpile groups 

Let r be a finite directed multigraph: it is defined by a finite set of vertices V and, 
for each pair of vertices v,w & V, a non-negative integer e{v,w), called the number 
of arrows from v to w. The total number of arrows going out of v is referred to as 
the outdegree of v and denoted by outdeg(f ); likewise, the indegree of v is the number 
of arrows going into v, denoted by indeg(v). (Computing these quantities, we do not 
take the loops, if any, into consideration.) We will assume that the graph F is strongly 
connected, that is, there is a directed path from any vertex v to any other vertex w. 
We will also suppose that our graph is Eulerian, that is, for every vertex v & V we 
have indeg(f) = outdeg(f). Under these assumptions, with the given graph T one 
can associate a certain finite abelian group >S'(r), called the sandpile group of F (see 
|Levt ILFj ). The group 5'(F) is defined uniquely up to isomorphism, and the simplest 
way to define it is through the Laplacian matrix of F. 

Let V = {f 1, . . . ,Vn}- The Laplacian matrix L = (lij) of size n x n is defined by 
its entries as 

la = - mdeg{vi), Uj = e{vi, Vj), if i ^ j. 

Let L C Z" be the lattice spanned by the rows of L. Evidently, L is a sublattice of 
Zq = {(«!, . . . , an) I «! + ■ ■ ■ + a„ = 0}. Then we set 

5(F) = Z^/L. 

It is known that the group 5'(F) can also be defined as follows. Delete any row 
and any column from matrix L and call the resulting (n — 1) x (n — 1) matrix L'. Let 
L' be the sublattice of Z"~^ spanned by the rows of L'. Then 5'(F) = Z"~^/L'. On 
a practical side, to compute the sandpile group, it is enough to reduce the Laplacian 
matrix by integral elementary operations on rows and columns to its Smith normal 
form, which is a diagonal matrix with integers {di, ...,dn-i,0) on the diagonal; then 
the group is ©"Ji^Zd^. 

In the next section we will define a series of Eulerian directed multigraphs F^ 
labelled by a prime number p and a natural number n and then show that, for p = 2, 
their sandpile groups are isomorphic to the automorphism groups of necklaces defined 
above. 



9 Generalized de Bruijn graphs 

Let F^ be the graph with vertex set F = Z„, the residues modulo n, and p directed 
edges from every vertex i to each oi pi, pi + 1, . . . , pi + p — 1. The outdegree of each 
vertex is thus equal to p. It is an easy exercise to check that the indegree of every 
vertex is also p and that the graph is strongly connected. Therefore, the sandpile 
group 5'(F^) is defined. 



We call these graphs generalized de Bruijn graphs, because the well-known de 
Bruijn graphs appear as a particular case V'^t,, see, e.g., jAdBj . 

The structure of the group ^(r^) for an arbitrary n is completely determined by 
the following two lemmas (Lemma [2] and Lemma [3]), the first of which treats the case 
of odd n and the second shows how to pass from any n to 2n. Before stating the 
lemmas, let us explain how one could actually arrive at the first, more difficult, one. 
Until the end of this section, we fix p = 2 and omit the superscripts 2 from various 
notations. 

Suppose that n is odd. Ideologically, the problem is quite simple: is suffices to 
find the Smith normal form of the integer matrix An explicitly defined by A„[0, 0] = 
An[n-l,n-l\ = -1, A„[0, 1] = An[n-l,n-2] = 1, An[i,i] = -2 ioT < i < n - 1 
and An[i, 2i] = An[i, 2i + 1] = 1 for z G Z„, i y^ 0, i ^ n — 1. The problem is purely 
technical, but rather difficult: to understand this, it is enough to look at the table 
of first 15 values of the sequence ^(rn) which coincides with the table of Section [6] 
and is quite non-trivial. The key difficulty is that, as a rule, there are three non-zero 
elements in each column and each row of this matrix, e.g. 



A, 



/-I 


1 




















0\ 





-2 


1 


1 























-2 





1 


1 




















-2 








1 


1 





1 











-2 











1 





1 


1 








-2 




















1 


1 





-2 























1 


1 


-2 





^0 




















1 


-V 



It would be much easier to treat a matrix where there are only two non-zero elements 
in each row and column. This goal is almost achieved through a trick invented by 
S.H.Chan in his Bachelor's paper |Ch] . 

Let us consider the operator given by the Laplace matrix of r„ in the basis Cq — ei. 



Cn— 1; (^n — ly 



that is, consider the matrix A' 



C„ ^ ■ A„- Cn, where C„ 



is a lower triangular matrix of size n with I's on the main diagonal and — I's on 
the adjacent diagonal. Multiplying a matrix by either C„ or C~^ is equivalent to 
elementary operations on its rows and columns, hence the Smith normal forms of 



matrices An and A'^ are the same. For the previous example we will obtain 



A' 



/o 























0\ 


1 


-2 


1 




















1 





-2 





1 














1 








-2 








1 








1 











-2 











1 





1 











-2 




















1 








-2 























1 





-2 





Vo 




















1 


-2/ 



We see that the lower-right minor of codimension 1 has the required property, and 
its Smith normal form can be found by drawing horizontal and vertical lines between 
the non-zero entries in each row and column of the matrix and considering the cycles 
obtained: 



-2 ■■■ 



























-2 



-2 



-2 



-2 



























1 



In this example we see two cycles of lengths 4 and 12 which are simply a visualization 
of the orbits of lengths 2 and 6 in the set Zg \ {0} under the doubling operator x i— )■ 2a; 
(in our case the orbits are {1, 2, 4, 8, 7, 5} and {3, 6}). It is readily verified that each 
orbit of length d adds a summand Z2d_i to the Smith group of such a matrix (where 
each row and each column contain one entry 1 and one entry —2), so for our example 
we obtain Zgs © Z3. Unfortunately, the presence of a nonzero first column spoils this 
clear picture, namely, it decreases the size of the group by a factor of n (more exactly, 
it leads to a subgroup of index n). For the example under study, any subgroup of 



index 9 is isomorphic to Z21. However, there are situations where such a group may 
have different subgroups of index n. S.H.Chan in a series of rather involved technical 
lemmas showed how exactly looks the resulting group S'(r„). In most cases, one must 
simply divide by n the order of the first cyclic group, corresponding to the orbit of 
the number 1. This is so for all odd integers up to 19. For n = 21, however, the set 
Z21 \{0} decomposes into five orbits {1, 2, 4, 8, 16, 11}, {3, 6, 12}, {5, 10, 20, 19, 17, 13}, 
{7, 14} and {9, 18, 15} of lengths 6, 3, 6, 2 and 3, respectively, but the sandpile group is 
actually equal to ZgQZ^, and not to ZlsZgQZ^ as one might infer from the previous 
rule. To state the exact formula proven by S.H.Chan, we need some notations. 

Let n be an odd number. For any element v G Z„ \ let l{v) be the length 
of its orbit under the doubling operator x 1— )■ 2x. Let if„ be the set of minimal 
representatives of all orbits. Now, denote by P„ the set of all prime divisors of n. For 
each g G P„ let q' stand for the maximal power of q that divides n and let q" = n/q'. 
Denote the set of all such residues q" by Vn- Finally, for an abelian group G and an 
integer k let kG be the subgroup of all elements of G of the form kx, x E G. Then 

Lemma 2. If n is odd, then the group ^(rn) has the following decomposition 



5(r„) - q'Z, 



2U'j")-l 



e 

veHn\Vn 



V2i(«)_l 



Example. Take n = 21. Then P„ = {3, 7}, q[ = 3, q'l = 7, /(7) = 2, q'^ = 7, q'^ = 3, 
/(3) = 3, H21 = {1,3,5,7,9}, V21 = {3,7}, so the first direct summand in the above 
formula is trivial, and the second gives: 



s{r 



21, 



^63 



z, 



'63 



Z7. 



The proof of the next lemma is much simpler: it follows from the fact that r2n is 
the directed line graph of r„ (see [Lev]). 



Lemma 3. Suppose that n = 2^m where m is odd. Then 

k 



S{Vn) = S{Tm) © Z 



m—1 

2fc 



^^2*+l-» 



i=2 



10 Circulant matrices 

An n X n matrix A = aij, i,j G Z„, over a field K is called circulant if its rows are 
the cyclic shifts of the first row, i.e. 

Ai+ij+i = Aij, foralH,jGZ„, 

where i + 1 and j + 1 are taken modulo n. 
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In particular, the permutation matrices associated with the powers of the cychc 
permutation (0, 1, . . . , n — 1) are circulant; they form a basis of the algebra C„(]K) = 
]K[Z„] of all n X n circulant matrices over K. We see that the algebra C„(K) has 
dimension n and is commutative. 

For a circulant matrix to be non-degenerate it is necessary (but not sufficient) that 
its rows (and columns) are aperiodic. Denote the group of non-degenerate circulants 
by C„(K)*. By the observation made above, it is commutative. In the case K = Fp, 
by studying the natural action of this group on the field Fpn considered as a vector 
space over Fp, it is easy to deduce that the Reutenauer group RG^ is isomorphic to 
the quotient C„(Fp)*/Z„, where Z„ is the group of permutation circulant matrices, 
those associated with the powers of (0, 1, . . . , n — 1). 

In the case p = 2, it was proved in |Ch] that the series of these quotient groups 
satisfies the same relations as those given in lemmas [2] and [31 The proof of these 
facts is not so involved as the proof of lemma |2] and relies basically on the primary 
decomposition theorem from linear algebra. The main theorem follows. 

Theorem. (S.W.Chan) For any natural n we have 

RGl - SiVl) . 

This theorem is quite remarkable, because it relates the objects coming from en- 
tirely different areas of mathematics. It is noteworthy that nobody knows any explicit 
isomorphism between the two groups in question, although the elements of both can 
be encoded by some sequences of O's and I's. 

In the paper |CHPj . this result is generalized to any prime number p as follows: 
RQP^ = S{TP) ©Zp_i. Moreover, that paper describes the structure of sandpile groups 
for the generalized de Bruijn graphs F^ for arbitrary values of p, not only prime. Of 
course, in the general case there is no analogue of the isomorphism theorem; however, 
if p is a power of a prime both groups are defined and the relation between them 
should be studied. 



11 Golomb's construction 

We conclude the paper with some experimental data related to another set of 1-1 
correspondences between the necklaces and irreducible polynomials mentioned above. 
Let a be a generator of the multiplicative group of the field Fg, q = p". S. Golomb 
|Gol] defined a bijection ipa '■ N^ -^ In which depends only on the orbit of the element 
a under the Galois group action, that is, tpa is completely determined by the primitive 
polynomial with one of the roots a. To a necklace z/ = (z/q, ..., i^^-i) we assign the 
following element of F^: 
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and then take its minimal polynomial, which we denote by ^a(z/). It is not hard to 
prove (see [Golj ) that the map ipa is 1-to-l for any a which is a root of a primitive 
polynomial, and that these maps are the same for all roots of one primitive polynomial, 
and are distinct for different primitive polynomials, so that they generate a group of 
automorphisms of necklaces whose order equals the number of primitive polynomials 
0(2*^ — l)/n, where is Euler's totient function. This group turns out to be abelian, 
too. 

Here is a table of these groups for p = 2 and 2 < n < 12 given together with the 
sizes of the orbits into which they split the set of necklaces (M„ stands for the order 
of the group): 



n 




M„ 


group 


orbits 


2 


1 


1 


Zi 


1 


3 


2 


2 


Z2 


2 


4 


3 


2 


Z2 


2+1 


5 


6 


6 


Ze 


6 


6 


9 


6 


Ze 


6 + 2 + 1 


7 


18 


18 


Z18 


18 


8 


30 


16 


Z2 0Z8 


16 + 8 + 4 + 2 


9 


56 


48 


Z2 © Z24 


48 + 8 


10 


99 


60 


Z2 © Z30 


60 + 30 + 6 + 2 + 1 


11 


186 


176 


Z2 © Zgg 


176 + 8 + 2 


12 


335 


144 


Z2 © Ze © Z12 


144 + 48 + 36 + 2 • 24 + 2 ■ 12 

+8 + 2-6 + 2-4 + 3-2 + 1 



Notice that the sequence of orders of these groups is not monotonous. As yet, 
nobody knows any relations between these groups and other mathematical objects, 
which was the case for the Reutenauer groups. 
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